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The induction of an s-wave component in a d-wave superconductor is considered. Near the {110}- 
oriented edges of such a sample, the induced s-wave order parameter together with d-wave component 
forms a complex combination d + e'^s, which breaks the time reversal symmetry (BTRS) of the 
pairing state. As a result, the spontaneous current is created. We numerically study the current 
distribution and the formation of the spontaneous flux induced by the current. We show that the 
spontaneous flux formed from a number of defect lines with {110} orientation has a measurable 
strength. This result may provide a unambiguous way to check the existence of BTRS pairing state 
at {110}-oriented boundaries. 



PACS numbers: 74.20.-z, 74.25.Jb, 74.50.+r 

A distinctive feature of a ci-wave superconductor is its 
sensitivity to inhomogeneities such as impurities, surfaces 
and interfaces. As one of the nodes of the d-wave en- 
ergy gap [Arf(k) ~ cos k x a — cosk y a in the momentum 
space] oriented along {110} direction is parallel to the 
normal vector of the surface, Andreev reflected quasi- 
particles from the surface may experience a sign change 
in the order parameter. Quantum constructive interfer- 
ence between the incident and Andreev reflected quasi- 
particle waves then leads to the formation of midgap 
states jjj, which can explain J2|,|| the zero-bias conduc- 
tance peak (ZBCP) observed when tunneling into the 
{110}-oriented thin films f|-|(|. Simultaneously, while 
the ci-wave order parameter is suppressed near the sur- 
face, a subdominant s-wave component of order param- 
eter may show up. It has been suggested that a broken- 
time-reversal-symmetry (BTRS) d + is pairing state |t]] 
could be generated near a surface PHlC}] or twin bound- 
ary ]il|-|l4j with {110} orientation. While the splitting 
of the ZBCP observed in a copper/insulator/Y-Ba-Cu- 
O (YBCO) tunneling experiment jl5) would be a signa- 
ture of BTRS pairing state fl(||l6| formed locally at sur- 
faces, other recent tunneling experiments have observed 
no such splitting ]T7|"[T9[] . Therefore, the existence of the 
surface BTRS pairing state in high-T c superconductors 
has not been fully established. Since it is a weak and 
secondary effect relative to the ZBCP in the tunneling 
experiment, one may argue that the samples pj|-^9[ not 
showing the splitting could be due to different surface and 
measurement conditions. Here we suggest to measure the 
spontaneous flux formed near the {110}-oriented edges 
as an alternative method to confirm the BTRS pairing 
state. In this paper, the spontaneous current distribu- 
tion or flow pattern in a square d-wave superconductor 
with {110}-oriented boundaries will be studied. From 
this information, we are able to calculate the sponta- 
neous flux distribution near the edges. We show that 
the total magnetic flux generated by a bunch of defect 
lines with {110} orientation should have a measurable 



strength. Its observation by experiments will unambigu- 
ously indicate the existence of the BTRS pairing state at 
the {110}-oriented boundaries. 

We first present a preliminary analysis based upon 
the Ginzburg-Landau (GL) theory to understand how 
a spontaneous supercurrent is generated by the surface 
BTRS pairing state in a <i-wave superconductor. The GL 
free energy for a two-component order parameter can be 
written as 

T = J d 2 r{a s \A s \ 2 + a d \A d \ 2 + /3 1 \A s \ 4 + 2 \A d \ 4 
+/33|A s | 2 |A d | 2 + / 3 4 (A: 2 A 2 + A 2 A* 2 )+ 75 .|nA:| 2 

+ 7d |nAS| 2 + 7 „[(n:A 5 n x A^ - n;A s n,A^) + c.c] 



where A* d are amplitudes of s- and d-wave components, 
o-s.d, Pi (i = 1,2,3,4), j s ,d,v ar e real coefficients, A is 
the vector potential, II = — iV — 2eA is the gauge in- 
variant momentum operator. Note that the free energy 
here is written in the coordinate system defined by the 
crystalline a and b axes (i.e., x and y axes in Fig. |l|) in 
YBCO. By varying the free energy with respect to the 
vector potential A, we obtain the current 

J = 2e{ 7s A s nA: + ld A d UA* 

+j v [A s U-A* d + A d n_A:]} + c.c , (2) 

where II = H x e x — Tl y e y with e XjV the unit vector along 
x (y) direction. We take the normal direction of the sur- 
face along the x axis and the parallel direction along the y 
axis. For a {100}-oriented surface where the x, y axes co- 
incide with the a, b crystalline axes, the two components 
of the order parameter are locked into a real combination 
d ± s by the mixed gradient term in the free energy. In 
this case, no spontaneous current exists. However, for 
a {110}-oriented surface where x and y are respectively 
along the y = x and y = —x directions as shown in Fig. [P, 
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this mixed gradient term vanishes and an imaginary com- 
bination d ± is is required to minimize the free energy. 
As a result, we find 

Jy = -4e 7 Jm(A s agA^ + A^A*) , (3) 

finite [J x — 4eIm(7 s A s 9xA* + ^dAddyA*^ is identically 
zero]. Away from the surface towards the bulk, the d- 
wave component saturates and the s-wave component de- 
cays to zero. Thus the parallel spontaneous supercurrent 
is limited to the surface region. The detailed value of 
Jy must be obtained by solving the GL equation for A* 
and A* d subject to the appropriate boundary conditions. 
However, the knowledge of Jy will not provide the infor- 
mation for the flow pattern since the supercurrent must 
circulate as loops in a realistic two dimensional system. 
The investigation of supercurrent flow is a challenging 
topic, it is difficult to employ the GL theory or other 
continuum theories to examine this problem. The ex- 
act diagonalization approach can attack this problem in 
an elegant way. In addition, the band structure and the 
short coherence length effects can also be incorporated. 
The weakness is that the finite size exact diagonaliza- 
tion method, if without modification, lacks the spectral 
resolution to resolve the resonances in energy. However, 
for the order parameter and the current (will be calcu- 
lated below) the finite size effect is small since these two 
physical quantities are determined by the local pairing 
interaction and local hopping integral, respectively. 

Within an extended Hubbard model, the Bogoliubov- 
de Gennes (BdG) [^0| equations can be written as: 

?(3 £)(?)-*(?)■ w 

where u" and u" are the Bogoliubov amplitudes at site i 
with eigenvalue E n , and 

#y = -%<5,j + V* - > ( 5 ) 

Au = Ao(i)*<j + A tf (i)<y i+ffj . (6) 

Here we choose the x and y axes to lie on the diagonals of 
the sample [as shown in Fig. [IJ], and 8 — ±x, ±y are the 
nearest-neighbor vectors along the x and y axes, t is the 
hopping integral, C7j is introduced to represent the impu- 
rity scattering potential. The energy gaps for on-site and 
nearest-neighbor pairing are determined self-consistently 

A (i) = V u?«r* tanh(E„/2T) , (7) 



A <5(i) = y E[<\h* + u r+K'1 ^HE n /2T) , (8) 

n 

where Vq and V\ are on-site and nearest-neighbor inter- 
action strength, respectively. Positive values of Vq and 



Vi mean attractive interactions and negative values mean 
repulsive interactions. In terms of the Bogoliubov ampli- 
tudes, the average current from site j to i is given by 

Jy = £{[/(£„KX + (1 - f(E n ))v?v?*} - c.c.} , 

n 

(9) 

where the prefactor comes from the spin degeneracy and 
f(E) = [exp(E/T) + l]" 1 is the Fermi distribution func- 
tion. 

We solve the BdG equations (0) self-consistently [ElJ . 
Throughout the work, we take V\ = 3t, fi = —t. This set 
of parameter values give A^ = 0.368£, T c — 0.569i, and 
the corresponding coherence length £ = hvp/irAd ~ 3a 
(In YBCO, £o ~ 15 A). As a model calculation, the value 
of hopping integral t will be adjusted to give a reasonable 
transition temperature. The amplitudes of the d- and the 
induced s-wave order parameters are defined as: 

A d (i) = ~[A*(i) + A_ 4 (i) - Ay(i) - A_y(i)] , (10a) 

A.(i) = J[A 4 (i) + + Ay(i) + A_y(i)] . (10b) 

In general, both components are complex, i.e., Ad, s = 
|A(j )S | exp(i<pd,s) with (f>d >a S [0, 2ir}. Their relative phase 
is defined as <f> = <j) s — 4>d- We first study the spatial 
variation of s-wave and d-wave order parameter and the 
relative phase in a square d-wave superconductor having 
{110}-oriented boundaries. The obtained current distri- 
bution is drawn in Fig. [I]. The current direction is de- 
noted by the arrows in the figure. This calculation is 
made under the open boundary condition for a clean sam- 
ple of size 10\/2a x lOv^a, T = 0.02i and V = -t. We 
find that the d-wave order parameter is suppressed near 
the boundary and approaches to the bulk value beyond 
a coherence length scale £o- The induced s-wave com- 
ponent near the boundary oscillates at an atomic scale 
and decays to zero in the bulk region at a distance of a 
few coherence lengths. Near the four corners of the sam- 
ple, the d-wave component is more strongly suppressed 
while the magnitude of induced s-wave component is en- 
hanced. Our numerical calculation shows that if the or- 
der parameter at site j within square OHAB (see Fig. |l|) 
is Ad — d r — idi and A s = s r + isi [for definiteness, d^j 
and s r ,i are taken to be positive], then at site j' within 
square OHGF which is related to site j by a ir/2 rotation 
around the origin O, the order parameter can be obtained 
as Ad = d r + idi and A s = — s r + isj. Correspondingly, 
the relative phase is 7r/2 + Scf) at site j and 7r/2 — Scf) at 
site j', where 8<j> = tan _1 (sj/s r ) — tan _1 (d r /dj) is site 
dependent. On the |x| = \y\ lines [the dotted lines in 
Fig. ^] , 8cf> — and the order parameter becomes exactly 
d + is. Therefore, the pairing state in a square d-wave su- 
perconductor with {110}-oriented boundaries is a BTRS 
d + e^s state with <j> varying spatially, and the existence 
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of the BTRS state is limited within a region of the or- 
der of £o near the sample edges. We have also found 
that the induced s-wave component of order parameter 
decreases with increasing temperature and on-site repul- 
sive interaction. For a given on-site interaction Vq = —t, 
the s-component or the surface pairing state vanishes at 
a temperature about 0.1T C [ pi) , hi our lattice model, 
the current flows into (or out of) each lattice site con- 
necting four bonds. The current conservation on each 
site is respected. As it can be seen from Fig. [I], the cur- 
rent flowing at four edges, AC, CE, EG, and GA, does 
not form a closed loop to surround the whole sample. 
Instead, the current distribution has a spatial reflection 
symmetry with respect to the x and y axis. Current flows 
pointing to the corners on the x axis but away from the 
corners on the y axis. At the left and right corners, the 
current flowing in the x direction is separated into two 
equal branches, which flow in the positive and negative y 
direction; while at the upper and lower corners, the cur- 
rent flowing in the y direction consists of two equal parts 
flowing in the positive and negative x direction. From an 
overall point of view, current flows clockwise or counter- 
clockwise in four triangles separated by the x and y axes 
[i.e., OAC, OAG, OGE, OEC]. Correspondingly, the total 
flux in two nearest-neighbor triangles for example OAC 
and OAG has opposite signs. In the upper-right inset of 
Fig. [j], we indicate the direction of flux in each triangle) . 
In contrast to the regular diamagnetic Meissner current 
which flows in a penetration depth, the spontaneous cur- 
rent is confined to a distance of a few £o from the sample 
edge because it is determined by the spatial variation of 
order parameter. 

To calculate the flux in each triangle, we use the 
Maxwell equation, VxB = (47r/c)J, which relates the 
magnetic field B to the current density J. By assum- 
ing that the lattice spacing between two consecutive su- 
perconducting layers is cq, the Maxwell equation can be 
written as 



where if ij is a horizontal link from left to right, <i>i and 
are respectively the magnetic flux passing through the 
upper and lower plaquettes with the link ij as a common 
boundary. The magnetic field and the corresponding flux 
is perpendicular to the lattice plane since both the cur- 
rent Jy and the vector connecting the centers of two pla- 
quettes are in the plane. The flux carried by a vortex is 
then obtained by summing over the flux encircled by each 
plaquctte in a triangle. Using the calculated current dis- 
tribution to solve Eq. (p"i|), the flux in the triangle AOC 
is found to be $ = £\ = 1.2£I , where C = 47ra 2 /cc 
and I = et/h. For t = 0.014 eV which gives T c ~ 92.4 K, 
and a — 4 A and c = 10 A, we have $ = 4.04 x 10~ 7 <& , 
where $o = 2.07 x 10~ 7 G • cm 2 is the superconducting 
flux quantum. Due to the small size we consider, the 



value is too small to be measurable. 

To measure the spontaneous flux as an evidence for 
the BTRS pairing state, we propose the following exper- 
iment. Create several long defect lines in a single crystal 
of d-wave superconductor. These lines lie parallel to the 
same {110} direction. To give a detailed analysis, we con- 
sider three parallel {110}-oriented defect lines in the sam- 
ple of size fOOv^a x 100v^a. The width of them is y/2a, 
\/2a, and 3-\/2a/2. The centers of the defect lines are 
located at 60a/V2, Wla/y/2, 159.5a/\/2. Therefore, the 
distance between the edges of two consecutive lines are, 
39a/ y/2, 56a/ \/2, respectively. In the numerical calcula- 
tion, we use the single-site potential of strength U = 100< 
to model the defect lines, and impose the periodic bound- 
ary condition to focus on the physical properties associ- 
ated with the defect lines. The spatial dependence of the 
order parameter and spontaneous current are shown in 
Fig. Ej. As expected, when the dominant d-wave order 
parameter is suppressed near the defect lines, the s-wave 
order parameter is admixed with a relative phase tt/2 
with respect to the d-wave order parameter. More im- 
portantly, the admixture of the s-wave and d-wave com- 
ponents is globally d + is: The relative phase is the same 
not only on the left and right sides of each defect line 
(i.e., d + is:d + is) but also near all defect lines. For the 
profile of the s-admixture near a defect line, a similar 
result that the relative phase is the same on both sides 
of a {110}-oriented boundary coupling two d-wave su- 
perconductors was obtained based on the quasi-classical 
theory [ |To| |. The same admixture of the s-wave compo- 
nent (i.e., d + is) on the two edges of a superconducting 
segment created by two consecutive defect lines, which 
comes from a small overlap of the electron wavefunctions 
at two defect lines, can gain more condensation energy. 
If the s-wave admixture became d + is at one edge and 
d — is at the other, there would be a 7r-phase jump in the 
middle region of the superconducting segment, which is 
physically unacceptable. More recently, it has also been 
found p2| that an s-wave component is admixed with the 
same sign on both sides of a {110}-oriented d-wave su- 
perconducting strip. This strip is in geometry similar to 
the superconducting segment we are considering here. As 
a consequence, the spontaneous currents flowing at the 
left and right sides of each defect line have opposite di- 
rections. In this respect, the system under consideration 
is different from an s-wave/normal-metal/d-wave super- 
conductor junction p3j| , where the BTRS pairing state 
appears in the normal metal, and the spontaneous su- 
percurrent has the same direction so that the generated 
magnetic flux has opposite sign across the corresponding 
region. In addition, our results show that the current 
profile near all defect lines is the same. With the above 
observation, the magnetic field induced by the current 
near these defect lines has the same direction. If the 
pairing state is d — is, the directions of the current and 
magnetic field are globally reversed. As shown in Fig. 0, 
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our results are independent of the width of and the dis- 
tance between the defect lines. Moreover, if the distance 
between the defect lines is of several coherence lengths, 
the contribution to the magnetic flux from each defect 
line is almost the same. With the chosen values of struc- 
ture parameter, the flux contribution from a single defect 
line of length 30 fj,m is estimated to be 3.4%$o- If twenty 
such defect lines are introduced into the sample, the total 
flux can be as large as 68%$o- In this sense, the effect of 
BTRS pairing state is greatly amplified whereas the in- 
formation from the tunneling spectroscopy is determined 
by the local property near one individual defect line or 
surface. We suggest to use the heavy-ion bombardment 
technique J24| to realize the desired arrangement of the 
defect lines. With this technique, the edge of the defect 
line can also be very sharp so that it has a well defined 
in-plane orientation. The magnetic flux can then be mea- 
sured by a SQUID microscope through a pickup loop |^5) . 
Here the pickup loop should be constructed to cover the 
central part of all defect lines. Since the proposed experi- 
ment does not involve as-grown grain boundaries or twin 
boundaries, the existence of a spontaneous flux would be 
a direct evidence for the BTRS pairing state. 
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FIG. 1. The current distribution in a square d- wave super- 
conductor with {110} boundaries. The current direction is 
represented with arrows. The value of parameters used in the 
calculation are T = 0.02t and Vo = —t. The upper-right inset 
shows the direction of flux in each triangle. 



FIG. 2. The spatial variation of the order parameter (a) 
and the spontaneous current (b) in a 100\/2a x 100\/2a sample 
containing three defect lines parallel to the {110} direction. 
The solid line (with filled square) corresponds to the d-wave 
component and the dotted line (with empty square) to the 
s-wave component. Here the distance is measured in units of 
a' — a/v2. 
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